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Abstract: Zelolitic imidazolate frameworks (ZIFs) are a new porous materials for CO; capture. In this work, 
we study the effect of H;O on CO; adsorption and CO;/N; separation in ZIF-8 and ZIF-90 using Grand 
Canonical Monte Carlo (GCMC) simulations. The results shows that the CO, uptake amount and CO;/N; 
selectivity in ZIF-90 are higher thant that in ZIF-8. In ZIF-8, the presence of H5O have no effect on CO; 
adsorption and CO;/N, separation. In ZIF-90, H2O has a cooperative effect on the adsorption of CO», the 
CO; uptake amount and CO»/N, selectivity increased with the presence of H5O. 


Key words: CO;; H20; Zelolitic imidazolate frameworks; Grand Canonical Monte Carlo simulations 


其 


收 稿 


HH: 2016-06-15 


起 金 项 
A 


第 一 作者 : 
Bh. d 


0 引言 


大 气 


: 国家 


K (199 


自然 科学 基金 资助 项 目 (No. 51676079); 湖北 省 
胡 建 ; 


d 


2-5, 9,1 


6 省 高 安 人 ， 博 RE 


然 科 学 基金 杰出 


! 青 年 项 目 No. 2015CFA046) 


主要 从 事 金属 有 机 骨架 吸附 分 离 CO» f 


文 ，liujing27@mail.hust.edu.cn 


P CO; 浓度 的 增加 所 导致 的 全 球 变 暖 是 人 
类 发 展 所 面临 的 一 项 重大 挑战 。 燃 煤 电厂 是 我 
大 的 Co; 排放 源 。 使 用 多 孔 吸 附 剂 从 电厂 


Hm 


气 中 捕 


^ 


3k CO,， 从 而 抑制 大 气 中 的 CO; 浓度 ， 被 认为 是 最 
为 有 效 的 CO; 减 排 方 法 之 一 贴 。 近 年 来 ， 多 种 多 孔 
材料 被 应 用 于 燃烧 后 CO» 捕获 研究 ， 其 中 金属 有 机 
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的 吸附 ， 而 Ni-MOF-74 的 CO 吸附 量 在 水 存在 的 情 
况 下 有 所 降低 。 然 而 ， 烟 气 中 的 HO 对 ZIFs 材料 
的 CO, 捕获 性 能 影响 的 研究 则 鲜 见 报道 中。 在 我 们 
前 期 研究 中 发 现 四 ，H20 的 存在 会 降低 ZIF-68 的 


CO; 的 吸附 量 ， 但 会 提高 材料 的 CON; 分 离 系数 。 


要 衡量 ZIFs 材料 在 燃烧 后 CO» 捕获 领域 的 应 用 潜 
力 ， 必 须 考察 烟 气 中 的 HO 对 ZIFs 材料 CO» 吸附 
性 能 和 CON» 分 离 性 能 的 影响 。 

本 文 针 对 两 种 具有 相同 拓扑 结构 但 包含 有 不 同 
有 机 链 的 ZIFs 材料 ，ZIF-8 和 ZIF-90， 采 用 巨 正则 
蒙特 卡 洛 (GCMC) 模 拟 的 方法 研究 烟 气 中 的 HzO 对 
这 两 种 ZIFs 材料 的 CO» 吸附 及 CON» 分 离 性 能 的 


影响 。 对 比 两 种 材料 在 COyN, 及 COJN;/H5O 两 种 


混合 气体 中 的 CO; 吸附 量 及 CON 分 离 系 数 ， 研 
究 H2O 对 ZIFs 材料 CO; 捕获 性 能 的 影响 ， 并 对 其 
影响 机 理 进 行 分 析 。 本 文 的 研究 结果 对 设计 满足 实 
际 工 况 应 用 条 件 的 ZIFs 材料 提供 了 理论 依据 。 


1 计算 模型 与 方法 
1.1 ZIFs 模型 

ZIF-8 及 ZIF-90 模型 中 各 原子 的 位 置 采用 XRD 
测定 的 实验 数据 区 ,9 确定 并 构建 模型 ， 模 型 如 图 1 所 
zn. ZIF-8 及 ZIF-90 具有 相同 的 SOD 拓扑 结构 ， 由 
金属 原子 Zn 与 不 同 的 有 机 链 链接 而 成 。ZIF-8 包含 
有 甲 基 咪 哗 有 机 链 ， 而 ZIF-90 则 由 醛 基 咪 唑 构成 。 
两 种 材料 的 结构 性 质 如 表 1 所 示 。 其 中 ， 可 接触 比 
表面 积 和 总 自由 体积 通过 Materials. Studio09 计 算得 
到 。 一 个 直径 等 于 No 动力 学 直径 (0.28 nm) 的 探 针 分 
子 用 于 计算 材料 的 可 接触 比 表 面积 (Sice)。 材 料 的 
总 自由 体积 为 未 被 骨架 原子 占据 的 体积 ， 因 此 直径 
为 0 ”nm 的 探 针 分 子 用 于 计算 材料 的 总 自由 体积 
(Vg). 
1.2 力 场 与 原子 电荷 

选用 Lennard-Jones 12-6 势能 模型 和 coulombic 
势能 模型 描述 ， 静 电 作 用 采用 Ewald 求 和 方法 描述 
气体 -气体 及 气体 -吸附 剂 之 间 的 相互 作用 力 : 
uy = Xie" -(3]« eus] 
式 中 : so = 8.8542 x 102 C2N-tm-2， 为 真空 电容 率 ; 
0op，Eo8 分 别 为 碰撞 直径 和 阱 深 。COs 和 N 分 子 为 
场 参 数 取 自 TraPPE 模型 L243。 在 该 模型 中 ，CO;, 为 
刚性 线 型 模型 ， 其 中 C-O 键 长 为 0.116 nm。 为 了 描 
R CO T BUE. C 原子 电荷 为 0.35e，O 原子 


键 长 为 0.11 
TEH, EWAN 原子 中 性 放置 一 个 带 有 +0.964e 


之 间 的 相互 作用 力 。 为 了 更 好 的 描述 气体 分 子 在 
ZIF-8 及 ZIF-90 中 的 吸附 行为 ，LJ 参数 取 自 优化 后 
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电 蓓 为 -0.7e。N, 分 子 为 三 点 刚性 模型 ， 其 中 N-N 
nm. N 原子 电荷 为 -0.482e， 为 平衡 分 


的 虚拟 原子 。H2O 则 取 自 TIP4P 模型 031， 该 模型 包 


括 有 三 个 电 蓓 原子 和 一 个 LJ 中心 原子 。 气 体 分 子 
力 场 参数 如 表 2 所 示 。 


表 1 ZIF-8 及 ZF-90 材料 的 结构 特性 
Table 1 Structural properties of ZIF-8 and ZIF-90 


Bd A 
ZIFs Y plgcm? | S«./m^?g  Viree/ cm?g 
数 /nm 
ZIF-8 1.699 0.924 2482 0.56 
ZIF-90 1.729 0.988 2267 0.52 


1 材料 结构 模型 (a) ZIF-8 (b) ZIF-90 
Fig. 1 Atomic structures of (a) ZIF-8 and (b) ZIF-90. 


ZIF-8 及 ZIF-90 力 场 参 数 取 自 UFF 72504, iz 
力 场 参 数 能 准确 描述 气体 分 子 与 ZIFs 材料 骨架 原子 


的 UFF 力 场 (9。 骨 架 原 子 力 场 参数 如 表 3 所 示 。 


ZIFs 材料 骨架 原子 电荷 取 自 通过 DFT 计算 拟 合 得 


到 的 ESP 电荷 03。 原 子 电 荷 如 图 2 所 示 。 


表 2 气体 分 子 L-J 势能 参数 及 原子 电荷 
Table 2 LJ Parameters and Charges for Guest Gas Molecules 


LJ 参数 
气体 分 子 c e/kb 电荷 
CO COC 2.8 27 0.70 
CO, O 3.05 79 -0.35 
N2 N N 33 36 -0.5075 
N5 com 0 0 1.015 
HO HMO H 0 0 0.52 


HO O 3.1536 78.03 0 
了 2O_com 0 0 -1.04 


表 3 ZIFs 材料 中 骨架 原子 工 -J 势能 参数 
Talbe 3 LJ Parameters for ZIF atoms 


原子 类 型 c £/kb 
Zn 2.338 43.084 
O 3.118 20.847 
C 3.259 36.483 
N 3.997 23.974 
H 2.440 15.288 


2 骨架 各 原子 电荷 (a) ZIF-8 (b) ZIF-90 
Fig. 2 Partial Atomic Charges of (a) ZIF-8 and (b) ZIF-90 


1.3 巨 正则 蒙特 卡 罗 模 拟 细节 

使 用 巨 正则 蒙特 卡 罗 (GCMC) 模 拟 方 法 计算 
ZIF-8 和 ZIF-90 材料 在 298 温度 条 件 下 的 单一 组 
分 气体 吸附 量 以 及 CON2y/H20 三 元 混合 气体 
组 分 气体 吸附 量 及 CO>/N;, 分 离 系 数 。 骨 架 中 各 原子 
为 刚性 ， 且 固定 在 各 自 位 置 。 考 虑 到 周围 环境 对 材 
料 吸 附 Co; 的 影响 ， 在 计算 中 使 用 三 维 周期 性 边界 
条 件 。 为 保证 计算 精度 ， 模 拟 盒子 选取 为 2x2x2 
个 单元 晶 胞 ， 截 断 半径 为 1.28 nm， 模 拟 步 数 为 2 x 
10“， 其 中 ， 前 1 x 10“ 用 于 平衡 体系 ， 后 1 x 107] 
于 热力 学 统计 研究 。 模 拟 计算 结果 为 绝对 吸附 量 ， 
而 实验 测量 结果 为 相对 吸附 量 ， 两 者 之 间 的 转换 公 
式 为 : Nex = Nabs 7 Vop, (2) 
RP, one 表示 相对 吸附 量 ，nsos 表示 绝对 吸附 量 ， 
所 表示 吸附 剂 孔 体积 ，ps 为 吸附 质 气 体 密 度 。pg 可 
通过 Peng-Robinson 方程 计算 得 到 。 和 气体 分 离 系 数 


L 


计算 公式 为 : 
S = (X1/X2z)(y2/y1) (3) 


Int 


中 x 和 y 分 别 是 相应 组 分 气体 在 吸附 相 以 及 气相 
的 摩尔 分 数 。 


N 


H 
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2 结果 与 讨论 


2.1 力 场 参数 的 验证 

在 对 模型 进行 分 子 模拟 研究 时 ， 选 择 合适 的 力 
场 参 数 对 于 模拟 结果 的 可 靠 性 有 决定 性 作用 。 本 文 
将 GCMC 模拟 所 得 到 的 ZIF-8 与 ZIF-90 的 CO; & 
NN; 吸 附 等 温 线 与 实验 所 得 到 的 吸附 等 温 线 相对 比 ， 
结果 如 图 3 所 示 。 结 果 表 明 模 拟 结果 与 实验 结果 很 
吻合 。 因 此 本 文 所 选取 的 力 场 参数 和 电荷 可 靠 。 


121 (a) ' J 
LoL 9 CO:ZIF-8 EXP il 
"ep — | —5— CO» ZIF-8 SIM 
© 0.8} ^ NoZIF-8 EXP ] 
E os L^ N: ZIF-8 SIM 
b: 
È 0.4 J 
EX E 
= 02 ] 
0.0 
0 20 40 60 100 
P/kPa 
2.5 T T T T 
) [9] 
29l o CO ZIF-90 EXP s | 


—e— CO» ZIF-90 SIM 
o CO»ZIF-8 EXP 
| -a— CO ZIF-8 SIM 


吸附 量 /mmol. g 
= 


0 20 40 60 80 100 
P/kPa 
3 模拟 与 实验 吸附 等 温 线 对 比 : (a) 298K, (b) 295K 
Fig. 3 Comparison of simulated and experiment adsorption 


isotherms: (a) 298K, (b) 295K 


2.2 单一 气体 吸附 

4 给 出 了 在 298 K lif CO,,. N% H50 fE 
ZIF-8 及 ZIF-90 上 的 吸附 等 温 线 。 其 中 CO, 压力 范 
围 为 : 0 kPa-100 kPa，N;, 压 力 范围 为 : 0 kPa-100 
kPa，H2O 的 压力 范围 为 0 kPa-10 kPa。 由 图 4 中 可 
以 看 出 ，CO, 在 ZIF-8 及 ZIF-90 中 的 吸附 量 明 显 高 
于 No XRH ZIF-8 及 ZIF-90 是 一 种 很 好 的 
CON; 分 离 材 料 。H2O 在 ZIF-8 中 的 吸附 量 很 低 ， 
这 是 因为 -CH; 官能 团 是 习 水 基 团 ， 会 阻碍 H2O 在 
ZIF-8 中 的 吸附 。 由 于 ZIF-90 中 的 -CHO 具有 亲 水 


(s 
a 
H 
[ep 
e 
BE 
N 
T 


性 ，H2O 在 ZIF-90 中 的 吸附 量 远 高 
90 中 的 吸附 量 。 


图 
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4 CO,, N: A H50 在 (a) ZIF-8 及 (b) ZIF-90 中 的 吸附 等 温 
线 
Fig. 4 Simulated CO2, Nz, and H,O adsorption isotherms in (a) 
ZIF-8 and (b) ZIF-90 


2.3 混合 气体 吸附 

为 了 研究 烟 气 中 的 HO 对 ZIF-8 及 ZIF-90 的 
CO 吸附 性 能 和 CON» 分 离 性 能 的 影响 ， 本 文 考虑 
了 两 种 不 同 的 气体 混合 物 : 1) COJN,. 2) 
COxNXYH2O， 为 模拟 实际 烟 气 ，CO; 分 压 为 15 
kPa, HO 为 5 kPa， 模拟 烟 气 总 压 为 100 kPa。 采 
用 N: 平衡 混合 气体 总 压 。 

图 5 给 出 了 ZIF-8 和 ZIF-90 在 两 种 混合 气体 
的 CO; 吸附 量 及 CON 分 离 系 数 。 从 图 5 中 可 以 
看 出 ， 在 COJN; 混合 气体 中 ，ZIF-90 的 CO; 吸附 
量 及 CO; 选择 性 要 高 于 ZIF-8 的 CO; 选择 性 。 这 是 
因为 ZIF-90 中 有 机 链 上 的 官能 团 -CHO 是 极 性 官能 
团 ， 而 在 ZIF-8 中 的 -CH3 官能 团 则 是 非 极 性 官能 
团 。 极 性 官能 团 相 比 非 极 性 官能 团 更 能 增强 CO» 在 
骨架 上 的 吸附 作用 能 09。 因 此 当 采 用 极 性 官能 团 取 
RET, CO 与 骨架 原子 之 间 会 具有 更 高 的 亲和力 ， 
从 而 会 增 大 材料 的 CO» 吸附 量 及 CON, 分 离 系 


L 


数 。 此 外 ， 我 们 还 可 以 发 现 HO 的 存在 对 ZIF-8 吸 
Iff CO, 及 分 离 COYN; 没有 影响 。 这 是 因为 H2O 在 
ZIF-8 中 的 吸附 量 极 低 ， 远 低 于 CO» 在 材料 内 的 吸 
附 。 另 一 方面 ， 尽 管 H2O 在 ZIF-90 中 的 吸附 量 相 
当 可 观 ， 但 HO 的 存在 却 对 材料 吸附 CO» 及 分 离 
COJN; 有 一 定 的 促进 作用 。 从 表 5 中 可 以 看 出 ， 
H20 的 存在 会 提高 CO, 在 ZIF-90 中 的 吸附 热 ， 因 
JE, HO 对 CO; 在 ZIF-90 中 的 吸附 具有 一 定 的 协 
同 吸附 效应 ， 从 而 提高 材料 的 COS 吸附 量 。 而 HO 
与 N? 之 间 则 缺乏 这 一 作用 ， 因 此 因此 H2O 的 存在 
增强 了 COYN2> 选择 性 系数 。 
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5 ZIF-8 和 ZIF-90 在 两 种 混合 气体 中 的 (a)CO, 吸附 量 和 
(b)CO;/N; 分 离 系数 


Fig. 5 (a) CO» uptake amount and (b) CO;/N; selectivity from 
two mixture in ZIF-8 and ZIF-90 


K 5 ZIF-8 和 ZIF-90 在 两 种 混合 气体 中 的 CO 吸附 热 
Table 5 CO; adsorption heat from two mixtures in ZIF-8 and 


ZIF-90 
ZIFs CO/N; CO;/N;/H;O 
ZIF-8 16.65 16.63 
ZIF-90 20.53 24.56 


本 文 使 用 巨 正则 蒙特 卡 罗 模 拟 (GCMC) 方 法 计 
算 了 HO 的 存在 对 ZIF-8 与 ZIF-90 的 CO; 吸附 性 
能 和 CON»; 分 离 性 能 的 影响 。 研 究 结 果 表 明 ，ZIF- 
90 的 CO; 吸附 及 COS/N» 分 离 性 能 要 强 于 ZIF-8。 在 
ZIF-8 "H, HO 的 吸附 量 远 低 于 COS 的 吸附 量 ， 而 
在 ZIF-90 中 ，H2O 的 吸附 量 高 于 CO» 的 吸附 量 。 
H20 的 存在 对 ZIF-8 的 CO; 吸附 性 能 和 CON» 分 离 
性 能 没有 影响 ， 而 在 ZIF-90 中 ，H20 与 CO, 之 间 
存在 协同 吸附 作用 。H20 的 存在 对 ZIF-90 的 CO; 
吸附 性 能 和 CON» 分 离 性 能 有 明显 的 促进 作用 。 
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